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Abstract 

In this paper, we present a statistical framework for modeling conditional quantiles of spatial processes 
assumed to be strongly mixing in space. We establish the Li consistency and the asymptotic normality 
of the kernel conditional quantile estimator in the case of random fields. We also define a nonparametric 
spatial predictor and illustrate the methodology used with some simulations. 
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1 Introduction 

Let {X, Y) be a pair of random variables with values in M'^ x IR and defined on a probability space (f2, A, P). 
Assume that the joint density of {X, Y) and the marginal density of X exist and are denoted respectively by 
f{x,y) and g{x). In the following, we suppose that F{-\x), the conditional distribution function of Y given 
X = X exists and we denote by f{-\x) the density of Y given X — x. For p e]0, 1[ and for fixed x G R'^, 
let fJ-p{x) be the conditional quantile of order p of F(-\x), that can be seen as a solution of the equation 
F(y\x) — p. Another alternative characterization of the pth conditional quantile (see for example Gannoun 
et al. [S]) is Hpix) = argmiuggR £; [(2p - 1){Y - 0) + \Y - 9\ \X = x]. 

We are interested to the non-parametric estimation of ^ip{x) in the case of spatial dependent observations. 
Nonparametric conditional quantile estimation technics have already been developed for non spatial (inde- 
pendent or mixing) real valued processes. Such results have provided useful tools for solving for example 
some prediction problems of strictly stationary processes satisfying the a-mixing condition. The existing 
results in the non-spatial case include the works of Matzner-L0ber [TU], Collomb H], Gannoun et al. [B], 
Laksaci et al. [S]. 

In nonparametric spatial estimation, the existing works concern mainly the estimation of a probability den- 
sity and regression functions, see the key references: Tran [T^], Biau and Cadre [5], Carbon et al. [3]. 
For the spatial quantile conditional estimation case, there exist only few results in our knowledge. Abdi et al. 
[T] considered the pointwise p— mean and almost complete consistencies of a double kernel quantile estimator 
for real-valued random fields. Hallin et al. [7] give a Bahadur representation and asymptotic normality re- 
sults of the local linear quantile estimator. Laksaci and Fouzia [5] consider the case where the regressor take 
their values in a semi-metric space and show the strong and weak consistency of the conditional quantile. 
In this paper, we will go beyond all these last spatial works and provide the Li consistency and an asymp- 
totic normality of a kernel conditional quantile estimate of a strictly stationary spatial process satisfying the 
a-mixing condition. In addition, we employe our results to solve some nonparametric prediction problems. 
The organization of this paper is as follows. The estimation procedure is presented in Section 2. Section 3 
gives some necessary conditions and then establishes the main asymptotic results. Section 4 is devoted to 
simulations results. Technical proofs are given in Section 5. 
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2 Nonparametric estimator of the conditional quantile 

Let us consider a strictly stationary process {{Xi,Yi),i £ X„) with values in M'' x M where {Xi,Yi) has the 
same distribution as (X, Y). For n e (N*)^, we define a rectangular region X„ by X„ = {i = (ii, . . . , i^v) G 
(N*)^, 1 < ife < rifc, fc = 1, . . . , iV}. We set n = ni . . . tin, and we write n — >■ cx) if minfc=i^...^jv nk oo. 
The well known kernel estimates of / and g are defined by 

/n(^' y) = J2 ^ i^ir^) ( Sr^ ) '''''^ •'^"^''^ " ^ ^ ( ) ' ^ '''''^ 

iex„ ^ ^ ^ ^ iein ^ ^ 

two probability density functions, and the bandwidths h = h(n) is a sequence of positive real numbers such 
that /i — > as n — > oo. The kernel estimate of the conditional density fn{y\x) is naturally defined by the 
ratio fn{x, y) over gn{x) while the estimator of the conditional distribution function (see the one introduced 
by Roussas [II]) is defined by 

T? t \ \ '^n{x,y) 



w 

— oo 



For a fixed x, the estimator of the pth conditional quantile noted fj.p.n{x) can be defined as the root of the 
equation Fn{z\x) — p. Alternatively, one can consider the local constant estimator defined by 



z.p,„(x) = argmin ^ {\Yi -6] + {2p - l){Yi - 9)) K 



X- Xi 



In this paper, we will focus on the study of ttie cisyniptotic belicivior of jJ^p^n- 

For the study of fp,n, one can 

adapt the technics developed in Zhou [14J. 



3 Main results 

To establish the asymptotic results, we will suppose that the sequence (Xi, Fi))ig(pj*-)jv satisfies the following 
mixing condition: there exists a function x '■ ^ with x{t) 4 as i — > cxd, such that whenever 
E,E' C (N*)^ with finite cardinals, 

a{B{E),B{E')) := sup {\P{A D B) ^ P{A)P{B)\; A e B{E), B e B{E')} 
< (t){CardE,CardE')xidist{E,E')), 

where B{E) (resp. B{E')) denotes the Borel tr-fields generated by {Xi,Yi)i^E (resp. {Xi,Yi)i^E'), Card E 
(resp. Card E') the cardinality of E (resp. E'), dist{E, E') the Euclidean distance between E and E' , and 
: — M"*" is a symmetric positive function which is non decreasing in each variable. Throughout this 
paper, we will assume that (/) satisfies 

(/) (n, m) < Cmin (n, m) , Vn, m G N (1) 



or 

0(n, m) < C(n + m + 1)" , Vn, m € N (2) 

for some k > 1 and some C > 0. If = 1, then the field (-'fi, yi)ig(N»)f< is called strongly mixing. In this 
paper, we consider the case where x(i) tends to zero at a polynomial rate, that is, 

Xi^)^0{^-P), (3) 
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with /? > 0. We fix a compact subset S of W^. Denote a = inf{y : F{y\x) > 0} and b = supjy : F{y\x) < 1}, 
we will suppose that V C [a,b] is a compact neighborhood of the unknown quantile For mixing 

coefficients with polynomial decreasing rate ([3]), the constraints on the bandwidth will be related to /3 by 
means of 

_ N{d+l){d + 2) + {d + l)l3 _ N{d + 2)-l3 
l3-N{d + 5) ' ^~ /3-N{d + 5)' 

_ N{d^ +M + 2) + {d+l)p _ N{d+l)^/3 

P-N{d + 4: + 2K,) ' '^^ l3-N{d + 4: + 2K)' 

Denote il„ = \/^f^- Let e be an arbitrary small positive number and set u{n) = J|^j^(logni)(loglogni)^+^. 
It is clear that X^nez" l/(nu(n)) < oo. In the sequel, we use the following hypotheses. 

(Al) / and g are respectively continuous on M'^+^ and M'^, g satisfies a Lipschitz condition, g{x) > 0, € S. 

(A2) There exists D > such that the pairs {Xi,Xj) and {{Xi,Yi), (Xj, Yj)) admit a density, say ^ij and 
/ij, as soon as dist{i,j) > D. Moreover, for some constant c > 0, 

\Aj{s,t)~ f{s)f{t)\<c, Vs,ieM'^+i 

and 

\gij{u,v)-g{u)g{v)\ < c, W,v£R''. 

(A3) i) /(*'^^(-, •) = g^tgyj (-, •) exists, is bounded and integrable for < i + j < 2, and < / < d. 
ii) F{y\x) has continuous second partial derivatives with respect to x. 

(A4) F{y\x) has continuous second derivative with respect to y. 

(A5) The kernel K is integrable, symmetric and is a lipschitzian density function on R"^ with compact 
support. Moreover J^^ \\s\\'^K{s)ds < oo. 

(A6) The kernel w is a symmetric and lipschitzian density function on M and has compact support. 

(A7) lim„^oo n/i''+2 (logfi)-^ = 0. 

(A8) The function fJ-p{x) satisfies a uniform uniqueness property on S: 

Ve > 0,3?? > 0,Vr : S R, sup |Aip(a;) - r{x)\ >e^ snp \F{^ip{x)\x) - F{r{x)\x)\ > rj. 

xgs xes 

(A9) nh^^{\ogny^ (w(n)) "-"f'^+'s) oo. 

(AlO) n/i''3(logn)''* (u(n)) oo. 

Comments on the hypotheses: 

Assumptions (A9) and (AlO) imply conditions (3.7) and (3.8) of Theorem 3.3 in Carbon et al. [3] and they 
also imply the classical condition nh'^'^^ / logn — oo. 

Assumption (A8) is introduced for getting consistency results on the quantile from those of the conditional 
distribution. 

In order to state the asymptotic results, we will suppose that (A9) and ^ or (AlO) and ^ are satisfied. 
The following two theorems give uniform almost sure convergence results of respectively F^{y\x) and fJ-p.n{x) 
and permit to establish the Li consistency of /Xp^n(x) (see Corollary [T]) . 
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Theorem 1 Assume (A1)-(A7) hold, then 



swpsvc[>\F^{y\x) - F{y\x)\ = O {Vt^) a.s. 



Theorem 2 // (A1)-(A8) are satisfied, the 



have 



sup |/ip,n(a;) - tJ-p{x)\ 0. 

Corollary 1 Assume (A1)-(A8) hold, then 

E [{Mp,n(Xn) - Mp(X„)}I{x„e5}] 0. 

To establish the following asymptotic normality of fip^ri{x) (Theorem|3]), we will suppose that for any {x, y) G 
S xV, there exists c > such that f(ii\x) > c. Moreover, we will assume that the following additional 
conditions on the bandwidth hold for some < 7 < 1. 

(CI) n/id(i+2W(i-7)) ^ 00. 

(C2) There exists a sequence of positive integers q — 00 with q = o (^(nft,'^*^^+^^(^~''')))^/^^) such that 

SESi ^"-'Xi^q) ^ and /j-'^^-^) »^-^(x(»))^-^ ^ 0. 

Theorem 3 Assume that (A1)-(A8), (CI) and (C2) hold. If there exists c > such that n/i'^+'* — >■ c, then 

in (x) - n (x)) 4 Af f^Eihlhl^^ a^(a:,Mp(x)) \ 
(Mp,n(x) A.,(a:)) ^JVyc ^(^^(^)i^) , if{f^,{x)\xW ) ' 



Bix,y) = 




92F(y|a;) 2 dg{x) dF{y\x) 



dxidxj 



g{x) dxi 



dxi 



\sfK{s)ds 



d^Fiy\x) 



t^w{t)dt 



a^{x,y) 



dy^ 

F{y\x) [1 - Fjylx)] 
9{x) 



(4) 
(5) 



3.1 Prediction 



Let (5i, i S Xji) be a K— valued strictly stationary random spatial process, assumed to be bounded, observable 
over a region In C and observed over a subset On of Xn, n = (ni, njv) S N^. The aim of this section 
is to predict ^j^, at a given fixed point io not in On C N^. In practice (e.g. for simplicity), we expect that 
depends only on the values of the process on a bounded neighborhood Vig C On- In other words, we 
expect that the process (^j) satisfies a Markov property, see for example Biau and Cadre |2J, Dabo-Niang 
and Yao [S]. Moreover, we assume that Vi^ = V + io, where V is a fixed bounded set of sites that does not 
contain 0. It is well known that the best predictor of ^i,-, given the data in Vi„ in the sense of mean-square 
error is 

i?(6ol6,ie ViJ. 

Let Vi = V + i = {u + i, u e V} for each i e N^, and d be the cardinal of V {d is also the cardinal of each 
Vi). To define a predictor of let us consider the R''-valued random variables = {^u, u e Vi C On}- 
The notation of the previous sections are used by setting Xi = Fi = ^i, i £ N^. 

As a predictor of fi^, we take the conditional quantile estimate = /Xp_n(Cio) order p, particularly the 
conditional median p — 0.5. We deduce from the previous consistency results, the following corollary that 
gives the convergence of the predictor f . 
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Corollary 2 i) Under the conditions of Corollary]^ we have 

E [{/ip,„(|io) - Aip(lio)}%i„es}] "-^^ 0- 
a) Under the conditions of Theorem\^ and ifnh'^^'^ — > 0, then 

where 

These consistency results permit to have an approximation of an 1 — a confidence interval of given 
by la = [a_(^iJ,a+(^io)], where 

a±(6o) = Mp.n(lio) ± Qi-^ ■ (6) 

ViIF/„(^p,„(Cio)|6o) 

where Q;; denotes the ^-quantile of the standard normal distribution, and the unknown parameters (of the 
asymptotic variance in Corollary [2]) are replaced by kernel estimates. 

Note also that the quantiles of order pi and p2 {pi < P2) can be used to construct a predictive interval that 
consists of the {p2 — pi)100% confidence interval with bounds /ipi,n(6o) Mp2,n(Cio)- 

4 A simulation study 

In this section, we study the performance of the conditional quantile predictor introduced in the previous 
section towards some simulations. Let us denoted by GRF{m, cr'^ , s) a Gaussian random field with mean m 
and covariance function defined by 



j2 



Set 

Xn = {i = {i,j) e (N*)2, 1 < i < 61, 1 < j < 61} (7) 
On = {i={i,j)e{N*f, l<i<2l, l<j<21}U{i = (22,j), l<j<15}. (8) 

We consider a random field (^i)igXn from the following model 

6 = C/i * (sin(2Xi) + 2exp{-(16Xi)2}) + Z;, i e (9) 

where X = {Xi)iex„ is a GRF{0, 5, 3), Z = (2'i)igx„ is a GRF{0, 0.1, 5) independent of X and U; = 
i X^jei "^^P (" ^^'2''^^ ) ■ choice of Ui in the model ^ is motivated by a reinforcement of the spatial local 
dependency. The field (^i, i G Xn) is observable over the rectangular region X„ and observed over the subset 
On defined in and dH]). 

We want to predict the values , . . . , ^i^^ at given fixed sites ii, . . . , not in ©„, with m = 10. The sample 
obtained from model observed in is plotted in Figure [T] below with the 10 non observable values of 
the field at ii, . . . , i„j. 
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Figure 1: The random field (^i) with non observed values ^i^, . . . in the white rectangular cases. 

As explained in Section |3.1[ for any k € {l,...,m}, we take the conditional quantile estimate ^i^. = 
/^pn(Cifc) a predictor of ^i^,, where ^i^. are observed on 0„ and the vicinity V = {—1,1} x {—1,1} or 
{-2,-l,l,2}x {-2,-1,1,2}. 

To compute /ip.n, we select the standard normal density as kernel K and the Epanechnikov kernel as w. For 
the bandwidth selection, we use the rule developed in Yu and Jones |13| . 

- /imean [^^^^T^ 

where /imean is the bandwidth for kernel smoothing estimation of the regression mean, </> and $, are respec- 
tively, the standard normal density and distribution function. 

To evaluate the performance of the predictor ^j^, we compute the mean absolute error (MAE): 

k=l 

The following Table gives the predictors of ^i,, k = 1,...,to forp e {0.05, 0.5, 0.95}, V = {-1,1} x {-1,1} 
on the left, V = {—2, —1, 1, 2} x {—2, —1, 1, 2} on the right and the prediction error. 




Table 1: Predictive data for V = {-1, 1} x {-1, 1} on the left and V = {-2,-1, 1, 2} x {-2, -1, 1, 2} on the 
right. 





p = 0.05 


True data 


p = 0.5 


p = 0.95 


p = 0.05 


True data 


p = 0.5 


p = 0.95 




0.1653 


0.2009 


0.1930 


0.2192 


0.1835 


0.2009 


0.2129 


0.2362 




-0.2553 


0.2315 


-0.2289 


-0.1862 


-0.2195 


-0.2315 


-0.1984 


-0.1766 




0.1516 


0.1966 


0.1990 


0.2362 


0.1912 


0.1966 


0.2129 


0.2362 




-0.5313 


0.4906 


-0.5062 


-0.4782 


-0.5472 


-0.4906 


-0.5313 


-0.5033 




0.2693 


0.2901 


0.2929 


0.3168 


0.2237 


0.2901 


0.2465 


0.2676 




-0.2748 


0.2535 


-0.2527 


-0.2289 


-0.2838 


-0.2535 


-0.2606 


-0.2401 




0.3696 


0.3941 


0.4007 


0.4269 


0.3805 


0.3941 


0.3834 


0.4269 




-0.5539 


0.5177 


-0.5295 


-0.5062 


-0.5472 


-0.5177 


-0.5313 


-0.5033 




-0.3678 


0.3217 


-0.3463 


-0.3193 


-0.3637 


-0.3217 


-0.3487 


-0.3231 




-0.2983 


0.2843 


-0.2702 


-0.2455 


-0.3096 


0.2843 


-0.2863 


-0.2671 


MAE 


0.0308 




0.0089 


0.0252 


0.0298 




0.0206 


0.0244 
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We derive from the results of Tables 1 a 90% predictive interval where the extremities are the 5% and 
95% quantiles estimates, for each of the 10 prediction sites (see Section [5TT] for more details). Note that 
these 90% predictive intervals contain the true values. The average length for the 10 intervals is 0.05604 for 
V = {-1, 1} X {-1, 1} and 0.04455 for V = {-2, -1, 1, 2} x {-2,-1, 1, 2}. 

The numerical results show that our proposed predictor gives good results on the above simulated field. A 
next step would be to apply the predictor to a spatial real data and deserves futur investigations. 

5 Appendix 

The letter C will be used to denote constants whose values are unimportant. Before proving the main results, 
let us give the following notations: 

/y ry 
f{x,z)dz, Ki{y) ^ w(z)dz. 
-oo J —oo 

Next, we need the following results. 

Lemma 1 Assume that (A3)-(A7) hold, then, 

supsup\Eipn{x,y) - 'ip{x,y)\ = O . 
yev xeS 

The proof of Lemma [1] is classical (see for example Matzner and L0ber jl0|). 

Lemma 2 Let Assumptions (Al)-(A6) hold. If (A9) and fip or (AlO) and (0) are satisfied, then, 

supsup|V'n(a;, y) - £^V'n(a;, y)| = O (^^n) as. 

Proof of Lemma [2j Define 

Ai{x,y) = Zi{x,y) - EZi{x,y), 

then Sn{x,y) = i)n{x,y) - Eipn{x,y) = ^ Ai{x,y). 

iex„ 

Lemma 3 Under Assumptions (Al), (A2), (A5) and (A6), we have 

E{Ai{x.y)f + J2 EAi{x,y)Aiix,y)^0 

iGX„ i#j 

The proof of Lemma |3] is similar to the proof of Lemma 2.2 in Tran [12]. 

Let us introduce a spatial block decomposition that has been used by Tran [12] and Carbon et al. [3]. 
Without loss of generality, assume that Ui — 2pqi ior 1 < i < N. The random variables Ai{x,y) can be 




7 



grouped into 2^qi x q2 x ■ ■ ■ x cubic blocks of side p. Denote 



(2jfc + l)p (2jfc + l)p 20„ + l)p 

Ai(x,2/), C/-(2,nJ,y)= ^ ^ Ai(x,y), 



fc=l TV 



ik=2jkP+l ijv = (2jN + l)p+l 

fc = l Af-l 



C/^(3,nJ,2/) = 

[/-(4,nJ,y) = 
and so on. Finally, note that 



(2jfc + l)p 2(ijv-l + l)p 

E E 



(2i«+l)p 

E ^i(^'y)' 



ik='^jkP+i ijv-i = (2jjv-l+l)p+l ijv=2iNP+l 

fc=l 



(2jfc + l)p 2(jjv-l + l)p 

E E 



2(iN+i)p 

E M^^y)^ 



ik=2jkP+l ijv-l = (2jjv-l+l)p+l ijv = (2iN + l)p+l 
fc=l,...,/V-2 



(2jfc+l)p (2iN + l)p (2ifc + l)p 

^ 5] Ai(x,2/), [/-(2^,nJ,y)= J] Ai(a;,y). 



ifc=2jfcp+l ijv=2jjvP+l 



jfc=2ifcP+l 
fc=l,...,iv 



For each integer 1 < i < 2 , let 



3k=0 

k = l,...,N 



Observe that, for any e > 

P{\Sn{x,y)\>e) 



E^"("'*'y) 



>e <2^P(|r-(n,z,y)| >e/2 



Without loss of generality, we consider just the case where i = 1 and we enumerate in an arbitrary way 
the q ^ qi . . . qpf terms U^{1, n, j, y) of the sum T^(n, 1, y) that we call Wi, . . . , W^. Note that U^{i, n, j, y) 
is measurable with respect to the cr-field generated by Zi, with i such that 2jkP + 1 < U- < {2jk + l)p, 
k=l,...,N. 

These sets of sites are separeted by a distance at least p and since the Zj are bounded, then we have for all 
z = |Wi| < Cinh'^y^p^WKWoo- Lemma 4.4 in Carbon et al. [3J ensures that there exist independent 

random variables , . . . , such that for alH = 1 , . . . , q, 

Em - W:\ < Cinh'')-'p''\\K\l^(^in,p'')x{p). 
Markov's inequality leads to 



pij2m- W:\ > £/2^+i < C2^+\nh'')-'p^q\\K\\^c^{n,p^)e~\{p). 



(10) 



\i=l 

By Bernstein's inequality, we have 
/ ? \ 



-eV(2^+i)2 



^ELiEW*^ + 2C{nhd)-^p^\\K\\^e/2N+\ 



(11) 
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Combining PH)) and ([TT|) . we have 

Pi\Sn{x,y)\>e) < 2^P(^|W^.-W,*| >£/2^+i I +2^P(|^W,*| >£/2^+i 



1=1 

2 



< 2--expi Z^W!!^ 



+C2^''+'^{n,p''){nh'')-'p''q\\K\\o.e-\{p). 
By LemmalU one has X]?=i EW*^ — O(g^) and since n = 2^p^g, we have 

P(|5„(a;,y)|>£) < 2^+^ exp | ^,N+,ct2^+^Cp^e } + ^^^^ V(a,P^)fe-'^£-^x(p). 

Let A > and set e = Afin, P — ^n^^^ ■ For the first part of Lemma [21 a simple computation shows that for 
sufficiently large n, 

■)N+l^^^f log" \ _^ noN+l^N,-cU-lc,-K 



P(|^„(x,y)|>A17„) < 2^^+Vxp^^5^^^5^^|^^^+C2^^+V^/i-'A-^17„^x(p) 



/3-2N 



< Cn-^ + C\-^h-'^nn'' , (12) 
with 6 > 0. Analogously, for the second part, as above we have 

P(|5„(x,y)|>Af^„) < 2^+^ exp | ^.m+^^cII^+^CX } + g2^^^n-/^-'^A-^»„^xb) 

< Cn-'' + CA-in''/i-''rJn" • (13) 

■^^^^ ] . Since 5 is a compact, it can be covered with cubes B^. 

having sides of length i?„ and center at Xk with c?„ < CR~'^. The compact set V can be covered with 1^ 
intervalls /; having length and center at yi, with ^„ < Cr~^ . We have 

dJn < C7ni/2/i-^ {h^+'n^y-' . (14) 
Define Sin = supj^gy^^axfe sup^gs^ \'>Pn{x,y) - i>nixk,y)\ 

S2n = maxmaxsup |V'n(a;fe,2/) - V'n(a;fc,yi)h 5'3„ = maxmax sup |£;Vn(a;fc, ?/;) - -EV'„(a;fc, j/)| 

k I y^J^ k I y^l^ 

Sin = sup max sup | £;?/;„ (a;fe, y) - E'ijjnix,y) \ , S^n = maxmax |7/;„(xfe, yj) - Eipn{xk,yi) \ ■ 

Then, we can write sup sup \ipnix, y) — Eipnix, y)\ < S'ln + 5'2n + S'sn + '5'4n + S^n- The proof of Lemma[2] 

yev xes 

follows easily from the combination of the two following lemmas. 
Lemma 4 Under Assumptions (A5) and (A6), 

S^n = 0{fln) a.s., i = 1,2,3,4. 

Lemma 5 Assume that (Al), (A2), (A5) and (A6) hold. If (A9) and (QP or A(10) and (0j are satisfied, 
then 

S511 = 0{fln) a.s. 
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Proof of Lemma [4l On one hand, since the kernel K satisfies the Lipschitz condition, we have clearly 

On the other hand, observe that h~^~^r-^ = (n/z^) = O (iln)- Since w satisfies the Lipschitz condition, 

which gives the proof of Lemma S) 
Proof of Lemma [5j For e > 0, we have 

P ^maxmax|V'n(a;fc,2/i) - -E?A„(a;fc,y;)| > < IndnP {\'ipn{xk,yi) - Eipn{xk,yi)\ > e) • 
Setting e = Ar2„ with A > and taking into account ([T^ and ([T5)l . it suffices to show that /nC?nn^''nu(n) — )• 

fH-2N 13- N 

and Indnh^'^^n " nu{n) 0, or Z„d„n''/i-^rj„" nu(n) ^ 0. 

First, observe that condition nh'^ oo implies that n > Ch^'^, so that 

g(d+l)(d+2)/2d (jf^-id+l){d+2}/2^ ^j^g-j 

Using dH]) and (fT5|) . we have 

;„d„S-''S7.(n) < C/i-('^+i)('^+2)/2(logS)-'^/2ni/2s-fH'*/2Hu(n) 
< Cn('''+3'^+i)/'^-''(logS)-'^/2M(n), 

which goes to if 6 > (d^ + 3d + 

Next, again (|14p and a computation show that 



/3~2iV _ 

Indnh'^n" nu(n) < C ri(logri)"(<i+=)-/3u(n)"(<i+5)-/3/i K(d+B)-/3 

which goes to by Assumption (A9) and (3 > N{d + 5). 
Analogously, (fT4|) and a computation show that 

^ JV(d+4+2n)-g 
2N -N{d^+id+2)-dl3) "' 

n(l0g n) «(d + 4+2K)-/3 u(n) N(d + 4 + 2K)-/3 /j JV(d+4+2»)-f) 



g-JV(d + 2) 2W -W(d+l)(d + 2)-dg) 



W(d + 5)-g 



which goes to by Assumption (AlO) and f3 > N{d + 4 + 2k). The conclusion of Lemma [5] follows from the 
Borel Cantelli's Lemma. 

Proof of Theorems [T] and [2j First, from Carbon et al. f3], we have sup3.g5 \gn{x) — g{x)\ = O (fin) a.s. 
Now, a standard decomposition gives 

sup sup \Faiy\x) - F{y\x)\ 

< — — - — -- I sup sup \iJnix,y) - i^{x,y) \ + sup sup (i^(y|a;)|5(a;) - ff„(x)|) 
ii^ixes 9n[x) ly£Vx£S yevxes 

^^ixGS gn[X) I y^VxGS y&V xeS xeS 
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Since by (Al), gn{x) is bounded away from 0, Theoreni[T]follows from the preceding inequality, Lemmas[T]and 
[5J Next, from (A8), to prove Theorem [21 it suffices to show that sup^.^^ \F{ijLp^n{x)\x) — F{iJ,p{x)\x)\ — > 0. 
We have 

\F{Hp,n{x)\x) - F{fJ.p{x)\x)\ < \F{flp^n{x)\x) - Fn{flp^n{x)\x)\ + |F„ (/Xp,n (a^) jx) - F{flp{x)\x)\ 

< \F{fip^n{x)\x) - Fn{fip,n{x)\x)\ < sup \Fniy\x) - F{y\x)\ . 

yev 

Thus, sup^gg \F{fip^n{x)\x) — F{^p{x)\x)\ < sup^gg sup^^^y \F„{y\x) — F{y\x)\ , so that Theorem [5] follows 
from an application of Theorem [1] 

Proof of Corollary [TJ First, by Lyapounov's inequality, we have 

E [^ip,n{Xn) ~ /ip(^„)] < E [|/ip,„(X„) - ^lp{Xn)\] < (E [{^ip^n{Xn) - Mp(^n))'])'^' , 

so, we can write 

[(Mp,n(^n) - Mp(X„)) I{x„es}] < {E[{tXp,^{Xn)-fipiX^)fl^x^^s}]y 
Then we also have 

£^[(Mp,n(^„)-A^p(^n))'l{X„eS}] < E 

An integration by parts gives 



.1/2 



sup {^ip,n{x) - fJ-pix)Y 

.xes 



< 2 



sup (/ip,„(a;) - ^J.pix)y 

xeS 



vP ( sup \fJ.p^ii{x) — lJ-p{x)\ > V j dv. 



Using Assumption (A8), we have that for n large enough, 

/ vP i swp\^p^nix) — ^ip{x)\ > v] dv < 2 / tP i snp\F{^p{x)\x) — Fn{^lp{x)\x)\ > t] dt 
Jo \x&S ' / Jo Vires' / 

< 2 [ tP (sup sup \F{y)\x) - F^{y\x)\ > t] dt, 
Jo \xesyev J 

and Corollary [1] follows from Theorem [TJ 
Proof of Theorem [Sj Since 

F(^[l.p(x)\x) =p = Fn{flp^n{x)\x), 

by Taylor's expansion, we have 



Mp,n(a;) - fipix) 



/n(/^;,„(a;)|a;) 



[F{np{x)\x) ~ Fn{np{x)\x)] , 



where fi*^(x) lies between fip^n^x) and ^ip{x). To prove the asymptotic normality, it suffices to show that 
the numerator is normally distributed, and the denominator converges to f {^,p{x)\x) in probability. We have 
the following propositions. 

Proposition 1 Under (Al)-(A6), (CI) and (C2), if there exists c > such that nh'^^'^ — > c, then 

Vm? {F^{y\x) ^ Fiy\x)) A M {cB{x,y),a\x,y)) , 
where B{x,y) and a^{x,y) are defined in ^ and ([5])- 
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Proposition 2 Assume that (Al), (A3)-(A7) hold. If (A9) and (QP or (AlO) and ^ are satisfied, then 



J ^J^^ sup \fn{y\x) - f{y\x)\ "4- 0. 

Proof of Proposition [H Proposition [1] is a consequence of the two following lemmas. 
Lemma 6 If Assumptions (Al) and (A3)-(A6) are satisfied, then 



EF^{y\x) - F{y\x) = —B{x, y) + o{h') + O 



where B{x,y) is defined in ([¥]). 

The proof of Lemma IH] is classical and therefore is omitted (see Matzner-L0ber |10j). 
Lemma 7 Assume that Assumptions (Al)-(A6), (CI) and (C2) hold, then 

Vm? (F„(y|x) - EFn{y\x)) 4 Af (O, a\x, y)) , 
where a'^{x,y) is defined in ([5]). 

Proof of Lemma [7l Assume for the moment that for any pair (ci, C2) G ffi^ with cf + ^ 0, 

Vfi^ [ci i9n{x) - Egnix)) + C2 (M^; y) - EiPnix, y))] 4 J\f (0, a*') , 
where cr*^ = [c1g{x) + c^il){x,y) + 2ciC2V'(a^j 2/)] / K'^{z)dz. Now, set 

1 



9nix) 

9{x) 



[■>Pn{x,y) - Eipn{x,y)] ^^"(^'^) ^^^(a;) _ Egn{x)] 

gn\X)Jlign[X) 

iPnix, y) - E^n{x, y)] - ^^^^^ [gn{x) ~ Eg^{x)] . 

9 [X) 



On one hand, according to (|16p . we have Vm?W* 4 Af (O, a^) , with being defined in ([5]). 
On the other hand, 



^[Wn-W*] 



'nl? 



1 



5n(a;) 9ix) 



[-4'nix,y) - Eil;nix,y)] 



■ip{x,y) E'ipn{x,y) 



9^ix) gn{x)Egn{x) 



(16) 



[gn{x) - Egn{x)] 



nh'^ 



which goes to in probability. 

Finally, we conclude with the decomposition 



Vn/i'i {F^{y\x) - EF^{y\x)) = Vn/i'^VF* + Vn/i'^ [W^ -W^] + 
Now, (|16p is proved following the same lines as the proof of Theorem 3.1 in Tran [12]. 
Proof of Proposition [2l We have 

sup|/„(y|x) - f{y\x)\ < sup \fn{x,y) - f{x,y)\ + sup f{y\x)\gn{x) ~ g{x)\. 

yev 9n[Xj yev 9n(X) yev 

p 

Hence to prove Proposition [21 it suffices to show that sup |/n(a;, y) — fix, y)\ — >■ 0, this is proved by following 
the same lines as in the proof of Theorem 3.3 in Carbon et al. [5]. 



12 



References 

[1] Abdi, S., Dabo-Niang, S., Diop, A. and Abdi, A. (2009). Consistency of a nonparametric conditional 
quantile estimator for random fields. To appear in Mathematical Methods of Statistics. 

[2] Biau, G. and Cadre, B. (2004). Nonparametric spatial regression, Statistical Inference for Stochastic 
Processes, vol. 7, pp. 327-349. 

[3] Carbon, M., Tran, L.T. and Wu, B. (1997). Kernel density estimation for random fields, Statistics and 
Probability Letters, vol. 36, pp. 115-125. 

[4] CoUomb, G. (1980). Estimation non parametrique de probabilites conditionnelles, C. R. Acad. Sci. Paris 
Ser I Math., vol. 291, pp. 427-430. 

[5] Dabo-Niang, S. and Yao, A. F. (2007). Spatial kernel regression estimation. Mathematical Methods of 
Statistics, vol. 16, Number 4, pp. 1-20. 

[6] Gannoun, A., Saracco, J. and Yu, K. (2003). Nonparametric prediction by conditional median and 
quantiles. Journal of Statistical Planning and Inference, vol. 117, pp. 207-223. 

[7] Hallin, M., Lu, Z. and Yu, K. (2009). Local linear spatial quantile regression, Bernouilli, vol. 15, 
Number 3, pp. 659-686. 

[8] Laksaci, A. Lemdani, M. and Ould-Said E. (2009). A generalized L^— approach for a kernel estimator 
of conditional quantile with functionnal regressors: Consistency and asymptotic normality. Statistics 
and Probability Letters, vol. 79, pp. 1065-1073. 

[9] Laksaci, A. and Fouzia, M. (2009). Estimation non parametrique de quantiles conditionnels pour des 
variables fonctionnelles spatialement dependantes, Comptes Rendus Mathematique, vol. 347, pp. 1075- 
1080. 

[10] Matzner-L0ber, E. (1997). Previson nonparametrique des processus stochastiques. These de I'Universite 
Montpellier II. 

[11] Roussas, G.G. (1969). Nonparametric estimation of the transition distribution function of a Markov 
process, Ann. Math. Statist, vol. 40, pp. 1386-1400. 

[12] Tran, L. T. (1990). Kernel density estimation on random fields. Journal of Multivariate Analysis, vol. 
34, pp. 37-53. 

[13] Yu, K. and Jones, M.C. (1998). Local linear quantile regression. J. Amer. Statist. Assoc. vol. 93, pp. 
228-237. 

[14] Zhou, Y. and Liang, H. (2000). Asymptotic Normality for Li Norm Kernel estimator of Conditional 
Median under a— Mixing Dependence, Journal of Multivariate Analysis, vol. 73, pp. 136-154. 



13 



